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Large pt processes at HERA, initiated by almost real and by virtual photons, provide information on the structure of the photon. 
We report on the latest measurements of dijets and large px particle production with the HI detector. This includes a leading order 
determination of an effective virtual photon parton density, of the gluon density of the photon, and comparisons with models. 



1 What is a photon? 

A photon is not just a photon and not just a hadron. In 
the leading order (LO) QCD framework it has two com- 
ponents, a direct one, which couples electromagnetically 
to the partons in the proton for example, and a resolved 
one. In its resolved state it fluctuates either into an on- 
shell qq-pair forming a vector meson (VDM part) or an 
off-shell qq-pair (anomalous part) , including any number 
of gluons, which may interact strongly with other par- 
tons around. In pictorial form this is expressed in Fig. |l|. 
In dijet or large pt particle production in ep collisions 
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Figure 1: 
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The dual nature of the photon and its associated nomen- 



we have two important scales to consider, the virtuality 
of the photon Q 2 and the mean Et of the hard jets or 
the large pr of an outgoing charged particle. In order 
to gain some intuition about the effect of the two scales, 
it is instructive to compare the life-times of the e7-state 
and the qq-fluctuation of the photon. Making use of the 
Heisenberg uncertainty relation, we find that the life-time 
of the e 7 -state falls like 1/Q 2 (i.e. from 0(10000) fm to 
O(100) fm for Q 2 from 0.1 to 100 GeV 2 . The life-time of 
the qq-fluctuation is constant for fixed pt such that the 
photon lives long enough do develop into a qq or even 
more complicated state as long as Q 2 < p\. In this situ- 
ation, in ep collisions, the partons in the proton are able 
to probe the partons in the photon target. However with 
increasing Q 2 , i.e. decreasing life-time of the photon, the 
photon becomes less resolvable. 

The hadronic structure of the photon has been well 
established both through measurements of the real pho- 
ton structure in 77 collisions at e + e~ colliders (PETRA, 
PEP, and LEP) and the measurements of jets in photo- 
production at HERA. Extending these measurements to 
the virtual photon structure is expected to provide new 
insight into the QCD framework, linking deep-inelastic 



(DIS), 7p, and 77 interactions. 

2 What is being measured at HERA? 

In LO QCD we expect contributions to dijet produc- 
tion from the direct processes of photon-gluon fusion and 
QCD-Compton and the resolved 2^2 parton processes 
shown in Fig. g. In ncxt-to- leading order (NLO) direct 
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Figure 2: LO direct (a) and resolved processes (b) 

and resolved processes can no longer be distinguished. 
Experimentally, the hard scale Q 2 and the fractional en- 
ergy of the photon, y, are determined from the measure- 
ment of the energy and angle of the scattered electron. 
The variable x^ ts , correlated with the fractional energy, 



of the parton in the photon participating in the hard 
scattering process is determined from the energy and di- 
rection of the two hardest jets and all of the hadrons in 
the event, using: 
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"Talk given on behalf of the HI Collab., ICHEP'98, Vancouver, 
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By cutting on x s ° ts , enriched event samples due to ei- 
ther direct or resolved processes can be obtained. At the 
parton level x^ ts = 1 for direct processes. 

The resolved ep cross section can be written as: 
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where the sum runs over all 2 — > 2 hard scattering matrix 
elements, folded by the density of parton j in the proton 
and parton i in the photon, and finally folded by the 
density of the photon in the electron. In the fractional 
energy, x p , of the parton in the proton, the measurement 
covers the range from about 0.01 to 0.1, where the parton 
densities of the proton are well known. 

The HI experiment has selected events in the range 
1.6 < Q 2 < 80 GeV| Land 0.1 < y < 0.7 in their most 
recent dijet analysis Uu. The jet selection was performed 
in the -f*p center of mass system (cms) using the inclusive 
k t algorithm!!! Events were required to have at least two 
jets satisfying the following criteria: 



hjoti - Jfietal < 1-0, 
E 2 > 30GeV 2 , and 



3.0 < f) < 

£ 1 T,jct2 
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< 0.25, 
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where fj and E t are the mean pseudorapidity (rj = 
— lntan(0/2)) and mean transverse energy of the two 
highest Et jets. The first two cuts make sure that the 
jets are confined to the acceptance of the detector, where 
they are well measured and that a; 7 is well determined. 
The restriction on the difference in r\ of the jets reduces 
the probability of misidcntifying a part of the photon 
or proton remnant as one of the high Et jets. The con- 
straints are such that the highest (second highest) Et jet 
has Et > 7 (4) GeV. This asymmetric jet selection will 
allow a comparison of the data with NLO calculations. 
With this selection HI obtains a sample of ~ 12000 dijet 
events for an integrated luminosity of 6pb _1 . 

The dijet cross sections measured as a function of 
x 7 , E 2 , and Q 2 have been corrected for detector accep- 
tance and resolution effects in separate ranges of Q 2 by 
applying an iterative Bayesian unfolding technique □ in 
x 1 and E 2 . The jet profiles and the pedestal energy out- 
side of thejets are reasonably well described Era by the 
HERWIG @ Monte Carlo £MC) used for the unfolding 
and also by the RAPGAPtl MC used to estimate uncer- 
tainties due to different models. A good description of 
the pedestal energy, also referred to as soft underlying 
event Era and in part caused by soft interactions of the 
photon and proton remnant, is of importance due to the 
steeply falling p^ distribution of the jets. 



For the virtual photon parton densities the Herwig 
MC uses the model by Drees and Godbole (DG) for 
suppressing the parton densities of real photons with 
increasing Q 2 by interpolating smoothly between the 
ln(p^/Ag CIJ ) behavior of real photons and the asymp- 
totic \n(p T /Q 2 ) dependence of the anomalous piece, i.e.: 



fq/j (-^7 J Pt> Q 

L(p T ,Q 2 ,t- 



2 ) = f qh (x 7 ,p 2 T ,0)L(p 2 T ,Q 2 ,u) 

2 ) = f gh (x J ,p T ,0)L 2 (p T ,Q 2 ,uj), where 
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where uj 2 is a free parameter to be chosen to describe the 
data. Another model of virtual photoa. parton densities 
is provided by Schuler and SjostrandQ, consisting of a 
vector meson dominance and a perturbative anomalous 
component with appropriate Q 2 evolution. 



3 Triple differential cross section 

The corrected triple differential cross section ia-shown as 
a function of x^ ts , E 2 , and Q 2 respectively intra. In each 
case the distributions are shown for ranges of the other 
two variables. Here, due to lack of space, we only show 
the cc J 7 ° ts distributions in Fig. |^. They can be seen to 
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Figure 3: The differential dijet cross section shown as a function of 
x^ tB for different regions of and Q 2 . Scale factors applied to the 
cross section are indicated. The error bars of the data points show 
the statistical and systematic errors added in quadrature. Also 
shown are model predictions (HERWIG) with 10% soft underlying 
event and two choices of the Q 2 suppression factor u). The direct 
processes as given by this model are indicated as shaded histogram. 
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peak towards x 3 ° ts = 1, where the direct contribution is 
expected to dominate. For a fixed Q 2 or life-time of the 
target photon the cross section shows a strong decrease 
with increasing E 2 of the probing jet. Increasing Q 2 , 
i.e. decreasing the life-time of the photon, while keep- 
ing the jet E 2 constant, leads to a decrease of the cross 
section and a diminishing relative contribution from re- 
solved photons (for x^ ts < 0.75). The HERWIG MC 
using GRV-LOcI for the parton densities of the real pho- 
ton and the DG model with to ~ 0.1 to 0.2 describes the 
data. In Fig. |I] the Q 2 dependence of the target photon 
is shown for a low (30 < E 2 < 45 GeV 2 ) and a high 
(65 < E 2 < 150 GeV 2 ) range in E 2 of the probing jet for 
a bin in £ J 7 ets dominated by resolved (0.3 < ir J 7 cts < 0.45) 
and by direct (0.75 < £ J 7 ots < 1.0) processes. We ob- 
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Figure 4: The differential dijet cross section shown as a function of 
Q 2 for a low and a high range in jet E 2 and for ranges in x^ ta dom- 
inated by either direct or resolved processes. Also shown are model 
predictions (RAPGAP) with a Q 2 suppression factor u> = 0.2. The 
sum of the LO direct and resolved contributions is indicated by the 
full line, the direct one by the dashed line. 



serve a strong decrease of the cross section with Q 2 . The 
resolved dominated contribution shows an even faster de- 
crease than the direct dominated one. This is what we 
qualitatively anticipate from the simple considerations in 
Sect. [|. They are supported quantitatively by the RAP- 
GAP MC as shown in Fig. |, which is in good agreement 
with the data. The stronger Q 2 suppression of the re- 
solved contribution is due to the additional Q 2 suppres- 
sion by the virtual photon parton densities. 

We conclude that the observed dependence of the di- 
jet cross section for x^ ts < 0.75 is consistent with that 
predicted for a resolved virtual photon with parton densi- 
ties evolving with Q 2 according to QCD motivated mod- 
els. 



4 Effective virtual photon parton density 

In order to determine the parton densities of the virtual 
photon, HI has adapted the single effective subprocess 
(SES) approximation, originally developed for use in pp 
collisionsEll and recently used to investigate real photon 
structure Ej. This approximation exploits the fact that 
the dominant contributions to the cross section comes 
from 2 — > 2 scattering matrix elements which have sim- 
ilar kinematic dependencies and differ mainly by their 
associated color factors. Therefore they can be replaced 
by an effective matrix element, M SES , and effective par- 
ton densities for the virtual photon, / 7 , and proton, f p . 
The resolved ep cross section of Eq. || then becomes 



SES f-f/e(UiQ ) fjix-yiPj^jQ 2 ) fp(x p ,p 2 -^) 
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where the effective parton densities are given by: 



N, 



A,p — yi(/g/7,p + fq/i,p) + Jal 
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and the sum runs over all quark flavors. 

For the extraction of the effective parton density only 
data with 0.2 < x^ ts < 0.7 and E^ jct > Q 2 , i.e. the con- 
dition for resolving the partons in the photon, are used 
for a second unfolding to correct the dijet cross section 
to the LO diparton cross section. This unfolding tries to 
correct for hadronisation effects and initial and final state 
QCD radiation. The systematic error Era includes those 
associated with the determination of the triple differen- 
tial cross section and additional errors arising in the sec- 
ond unfolding and amounts to an average error of ~ 40%. 

The resulting effective parton densities-divided by 
the fine-structure constant a, are shown intra as a func- 
tion of £ 7 , pip, and Q 2 respectively. In x 1 they show a 
small rise towards high they are consistent with 

a scaling behavior ~ Inp 2 -., and in Q 2 they are suppressed 
with increasing virtuality of the photon as predicted by 
QCD. They are reasonably well described by GRV-LO 
and the Drees and Godbole model with u ~ 0.1 and the 
parametrisations SAS-1D and SAS-2D of the model by 
Schuler and Sjostrand. In Fig. ^| we show only the re- 
sults as a function of Q 2 for a fixed scale p\ = 85 GeV 2 
of the probing jet and two different bins in x~ i . Also 
shown are two photoproduction data points from HI EH. 
The Q 2 evolution of the data are compared to the two 
models mentioned above and to a single p-pole suppres- 
sion factor, (m 2 /(m 2 + Q 2 )) 2 , characteristic of a simple 
VMD model. The p-pole factor clearly underestimates 
the data, while the logarithmic suppression, as expected 
from QCD with decreasing life-time of the photon, is in 
agreement with the observation. 



3 



H1 PRELIMINARY 



— Rho pole 

— GRV-DG(oi = 0.1) 




Q 2 /GeV 2 

Figure 5: The LO effective parton density of the photon, divided 
by the fine-structure constant a, as a function of photon virtuality 
Q 2 for p^, = 85 GeV 2 and two values for x 7 . The error bars of the 
data points indicate the statistical and systematic errors added in 
quadrature. Also shown are predictions from the DG model with 
u) = 0.1, using GRV-LO parton densities for real photons and the 
SAS-1D (dashed line) and SAS-2D (dot-dashed line) parametrisa- 
tions. The steeply falling solid curve shows the photoproduction 
data point extrapolated by a p-pole factor. 



5 Gluon density of the photon 

Measurements of the photon structure in e + e _ interac- 
tions are sensitive to the quark structure of the photon 
and only indirectly through the QCD evolution to the 
gluon structure. The data so far have not been precise 
enough to allow a determination of the gluon density of 
the photon. Recently, studies of photo-produced dijets 
and high px charged particles at HERA have shown that 
these data are sensitive to both the quark and gluon con- 
tent of the photon. A first measurement by HI of the LO 
gluon density E3 showed that it is not large at high x 7 , 
in contrast to the LAC3 t3 parametrisation, and does 
not have a_steep rise towards low x 7 as predicted by 
the LAC1E3 parametrisation. The minimum x^ ts value 
which was reached in this measurement was 0.04. It is of 
interest to reach x 1 values as low as possible, since most 
suggested gluon densities are predicted to rise. Experi- 
mentally it is rather difficult to reach low x 7 . Because of 
the relation x 1 ~ E T exp" 7 * /{2E 1 ) for Q 2 ~ GeV 2 , low 
x 1 implies small E^ (in conflict with a good correlation 
with the hard dipartons), large rj (soft underlying event 
and detector acceptance) , and large E 1 (decreasing event 
rate). 



Hl_.has contributed two different analyses, one on 
dijets Ej Jisihich is still in progress axid one on high px 
particles E3 which has become final E3 after the confer- 
ence. The latter analysis does not suffer from the energy 
scale uncertainty of the calorimeter and is less sensitive 
to the soft multiple interactions compared to the dijet 
analysis, but has the drawback of stronger sensitivity to 
uncertainties in the fragmentation. 

Both analyses require the scattered electron to be 
tagged, which restricts the photon virtuality to Q 2 < 
0.01 GeV 2 . The energy fraction y of the radiated pho- 
ton is required to be in the range 0.5 (0.3 for the high px 
analysis) < y < 0.7. In the dijet analysis at-least two jets 
have to be found using a cone algorithm E-3 with a cone 
radius of 0.7. Additional requirements on the jets are: 
prjet > 4 GeV, M jen , jet 2 > 12 GeV, -0.5 < r} iet < 2.5 
and |?7jcti — ^7jct2 < 1 in the HERA-system. Comparisons 
of data and MC can be found inli-i The corrected dijet 
cross section as a function of a; J 7 ts , after unfolding □, is 
shown in Fig. ^. It is compared to LO predictions by the 
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Figure 6: The ep dijet cross section for Q 2 < 0.01 GeV 2 as a 
function of x^ ts . The total error bars on the data points reflect 
the statistical and systematic errors added in quadrature. The 
data are compared to a LO model (PHOJET) using the GRV-LO 
parton densities for the photon and the proton. The contributions 
from direct and resolved, quark or gluon initiated, processes are 
indicated. 

PHOJET0 MC using the GRV-LO parton densities for 
the photon and the proton. The transverse energy £x 
of the jets is used for both the renormalisation and fac- 
torisation scale. Also indicated are the different contri- 
butions from the direct and resolved processes; the latter 
are split into contributions initiated from either a quark 
or a gluon on the photon side. It is clear from the fig- 
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lire that there is a sizeable contribution from gluons and 
that the data are precise enough to constrain the gluon 
density with a precision of ss 30%. 

In the high pt charged particle analysis the further 
requirements are: tracks which have a pt > 2 GeV and 
1 77 1 < 1 in the HERA laboratory frame. In Fig. 0a the 
corrected pr distribution and fits to the HI and pp data 
are shown. A QCD inspired power-law expression of the 
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Figure 7: a) The inclusive -yp cross section for charged particles 
as a function of px ■ The error bars on the data points denote the 
statistical and systematic errors added in quadrature. The curves 
indicate power-law fits to the HI, CDF, and UA1 data; the latter 
have been normalised to the HI data at px = 2 GeV. b) The ratio 
of data to an NLO QCD calculation with scale §p T - 



form A{l+p T /p Ta )- n was fit. The fit gives A = 5.44 ± 
0.66 mb and n = 7.03 ± 0.21 (stat.+syst.) and describes 
the data well over the whole range in pt- Similar fits to 
pp data from UAlllj and CDFcj have been normalised 
to the 7p cross section at pt = 2 GeV. The high pt tail 
in the jp data is clearly larger than in pp collisions at 
similar cms energies. This can be understood as being 
due to extra contributions in jp, namely the direct and 
the point-like resolved component. Similar conclusions 
have been obtained for 77 collisions by OPALcJ. 



The ratio of data to theory, as given by an NLO cal- 
culation including direct and resolved contributions E3, 
and using as scales ^ = ^ = ^ = £ p 2 h ^ with 
£ = 1 /2, 1, and 2, is shown in Fig. 0b. The pseudorapidity 
distributions for data and the NLO calculation for pt > 2 
and pt > 3 GeV are given in Fig. ||. The data are well de- 
scribed by NLO, in particular for £ close to 1. However, 
one finds that the NLO prediction is rather sensitive to 
the choice of scale (£). This effect is smaller for larger pt- 
For each event with at least one charged particle with 
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Figure 8: The inclusive 7p cross section as a function of 7/ for 
charged particles with a) px > 2 GeV and b) px > 3 GeV in 
comparison with an NLO QCD calculation with scale £p T - The 
error bars indicate the statistical and systematic errors added in 
quadrature. 
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Pt > 2.6 GeV, the variable x^ cc = J2pt ex P( — v)/E- 
where the sum runs over all tracks with itr > 2 GeV, 
is calculated. It shows a good correlationEijto the true 
Xj. To obtain x 7 an unfolding procedure E3 was used. 
Then the LO gluon density was unfolded allowing for 
an uncertainty in the quark densitiespof ~ 15% as given 
by three different parametrisationsQ'EHi3 and using differ- 
ent hadronisation models. The resulting gluon density 
(only available after the conference) is shown in Fig. || 
as a function_of x 7 and compared to an older extraction 
using dijets t3 and to three different parametrisations. 
The mean p T of the hard scattering process for this data 
sample is 38 GeV 2 according to MC, which was used as 
the scale for the comparison with three narton density 
parametrisations, GRV-LO □, SAS-ID-LOB, and LAC1- 
LOLL3. The mean p\ for the dijet sample is 75 GeV 2 . The 
results confirm that the contribution of the gluon to the 
photon structure is significant. The gluon density rises 
with decreasing x 1 and is best described by GRV-LO. 
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Figure 9: The LO gluon density in the photon from charged high 
px tracks (full circles) and from dijets (open circles) with a mean 
scale p^, = 38 GeV 2 and p^, = 75 GeV 2 respectively. The full error 
bars give the statistical and systematic errors added in quadrature. 
The lines show the predictions of three different parametrisations. 



6 Summary 

At HERA, with dijet and high px particle production, 
we can probe the virtual photon structure over a large 
range in the scale ~ p T and in the life-time of the photon 
~ 1/Q 2 . 

In the measured dijet cross section as well as the 
extracted effective virtual photon parton density, we ob- 
serve the expected logarithmic suppression with Q 2 , as 
the life-time of the photon decreases, of the resolved pho- 
ton contribution. 

Using two complementary methods, high px charged 
particles and dijets in almost real photoproduction, the 
leading order gluon density of the photon was determined 
and found to be rising with decreasing x 1 . 

So far we have only scratched the surface in that we 
have found consistency with rather global expectations 
and in that we have learned how to do leading order 
determinations of parton densities. There is a lot more 
to do, stay tuned for 1999. 
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